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ABSTRACT 

Context. The prompt emissions of gamma-ray bursts (GRBs) are seeded by radiating ultrarelativistic electrons. Kinetic energy domi- 
nated internal shocks propagating through a jet launched by a stellar implosion, are expected to dually amplify the magnetic field and 
accelerate electrons. 

Aims. We explore the effects of density asymmetry and of a quasi-parallel magnetic field on the collision of two plasma clouds. 
Methods. A two-dimensional relativistic particle-in-cell (PIC) simulation models the collision with 0.9c of two plasma clouds, in the 
presence of a quasi-parallel magnetic field. The cloud density ratio is 10. The densities of ions and electrons and the temperature of 
131 keV are equal in each cloud, and the mass ratio is 250. The peak Lorentz factor of the electrons is determined, along with the 
orientation and the strength of the magnetic field at the cloud collision boundary. 

Results. The magnetic field component orthogonal to the initial plasma flow direction is amplified to values that exceed those ex- 
pected from the shock compression by over an order of magnitude. The forming shock is quasi-perpendicular due to this amplification, 
caused by a current sheet which develops in response to the differing deflection of the upstream electrons and ions incident on the 
magnetised shock transition layer. The electron deflection implies a charge separation of the upstream electrons and ions; the resulting 
electric field drags the electrons through the magnetic field, whereupon they acquire a relativistic mass comparable to that of the ions. 
We demonstrate how a magnetic field structure resembling the cross section of a flux tube grows self-consistently in the current sheet 
of the shock transition layer. Plasma filamentation develops behind the shock front, as well as signatures of orthogonal magnetic field 
striping, indicative of the filamentation instability. These magnetic fields convect away from the shock boundary and their energy 
density exceeds by far the thermal pressure of the plasma. Localized magnetic bubbles form. 

Conclusions. Energy equipartition between the ion, electron and magnetic energy is obtained at the shock transition layer. The 
electronic radiation can provide a seed photon population that can be energized by secondary processes (e.g. inverse Compton). 

Key words, gamma rays: bursts - acceleration of particles - shock waves - magnetic fields - ISM: jets and outflows - meth- 
ods numerical 



1. Introduction 

1.1. Observations and Context 

Gamma ray bursts (GRBs) are eruptions of electromagnetic radi- 
ation at cosmological distances. One group of GRBs, those with 
a long duration, is attributed to the implosion of supermassive 
stars. This is suppo rted by ob servations, where GRBs precede 
supernovae (Hjort het al.|2003| ) and of particularly viol ent stellar 
explosions that show some resemblances with GRBs (Kulkarni 
|et al.|1998 ). GRBs are thought to be signatures of plasma ejec- 
tion from a forming compact object, such as a neutron star or a 
black hole. 

The fireball mod el due to |Meszaros & Rees] ([1992) and |Rees| 
& Meszaros (1994]) assumes that the plasma is ejected in form 
of a highly relativistic collimated jet by extreme supernovae (hy- 
pernovae) and that the jet dynamics is kinematically driven. It 
has been used to explain the anisotropic radiation bursts. In this 
model, plasma clouds collide due to the nonuniform flow speed 
and density of the jet, which is a consequence of the nonsta- 



tionarity of the jet source. The clouds can collide within the jet 
with Lorentz factors of a few and the cloud densities can proba- 
bly vary by about an order of magnitude. The resulting internal 
shocks move with Lorentz factors of a few through the jet. They 
are thoug ht to be responsible for the observ ed prompt phase of 
[T9991 |Fox & Meszaros|[2(306| ). The prompt emis- 
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sions due to the jet thermalization precede the longlasting after- 
glow, which has its origin in the interaction of the jet plasma 
with the ambient medium. 

The underlying mechanisms causing the observed electro- 
magnetic radiation are still not fully understood. Highly polar- 
ized gamma-ray emission suggests th e presence of an ordered 
magnetic field in the emitting zone ( [Steele et~a l. 2009). These 
primordial magnetic fields, which can be amplified further by 
the internal shocks, together with ultrarelativistic electrons give 
rise to electromagnetic emissions. The resulting photon seed can 
be upscattered to higher energies by secondary processes ( Kirk 
& Re ville 2010). It is not yet clear if, and to what extent, sta- 
ble and large-scale magnetic fields are generated or amplified by 



2 



G. C. Murphy et al.: Field amplification and electron acceleration in a protoshock 



instabilities close to internal sh ocks ([Medvede v & Loeb||1999| 
Brainerd 2000; Waxman 2006). The presence of ultrarelativis- 



tic electrons in the jet also cannot be taken for granted. The 
dominant blackbody radiation component of GRB jets suggests 
a plasma temperature of -100 keV (Ryde 2005 ). The majority of 
jet electrons are thus only moderately relativistic. The energetic 
electrons responsible for the nonthermal radiation component of 
the prompt emissions (Ryde 2005) must consequently be accel- 
erated within the jet, probably by the internal shocks. 

The likely involvement of plasma collisions for both mag- 
netic field generation/amplification and for electron acceleration 
and the possible presence of a guiding magnetic field, motivates 
our simulation study of the early stages of plasma cloud collision 
and shock formation. 



1.2. General Plasma behaviour 

The single- fluid magnetohydrodynamic (MHD) approximation 
can be used to model the ejection of relativistic jets by compact 
objects and to examine their time evolution (Nishik awa et al.| 
2005). However, it can not adequately describe the plasma dy- 
namics within the shock transition layer, in which the ultrarel- 
ativistic electrons and the strong magnetic fields are generated. 
Many different wave-modes (e.g. upper hybrid waves, whistler 
waves and Bernstein mode waves) and instabilities (Bre t|2009] ) 
not captured by a MHD model are present and they interplay. 
The wave and instability spectrum depends critically on multi- 
ple parameters of the bulk plasma, among others the electric and 
magnetic field orientation and strength, the plasma composition, 
the pre- and post shock density ratio and the temperature. In or- 
der to model the acceleration of particles, the generation of mag- 
netic fields by internal shocks in GRB jets and their small scale 
structure, it is necessary to use a kinetic particle-in-cell (PIC) 
simulation code. 



1.3. Modelling and PIC simulations 

The theory of collisionless magnetised shocks divides natu- 
rally into shocks with a quasi-parallel magnetic field (treated in 
this work) and those with a quasi-perpendicular shock. Quasi- 
perpendicular shocks have been well studied in the past and in 
the context of SNR or Solar system shocks using both analyt- 
ical and numerical approaches. Single-fluid analytical models 
have been unable to describe successfully the TeV emission from 
SNRs but hybrid or PIC simulations may provide better insight 
( Kirk & Dendy|200"T] ). Nonrelativistic numerical models of per- 
pendicular and parallel electrostatic shocks have been devised 
using hybrid methods, for example by ( [Leroy et al.| 1981 [[Q uest 



1988]). Perpendicular shocks have been modelled by|Lee et al. 
( 200 4);|Scholer & Matsuki^(|200^;|Amano & Hoshino| ( |2007| ); 



Umed a et al.| ( |2009| ); |Lembege et al. (2009) and many more au- 
thors with PIC simulations, while |Sorasio et al . ('2006') addressed 
fast unmagnetized shocks. Oblique, strongly magnetised shocks 



have been stud i ed by Lembeg e & Dawson] (|1989|); |Bessho & 
|Qhsawa|(^999l);|Dieclmiann et al.| (|20Q8|); |Sironi & Spitkovsky 
( |2009| ); |Shikii & Toida| ( |2010| ); |Murphy et al.| ( |2010a| )~ 



Considerable work has been done modelling interpenetrat- 
ing and colliding plasma streams in the context of GRB jets, 
with t heir a ttendant wave modes and instabilities. Frederiksen 
et al. ( 2004| ) studied using 3D PIC simulations two initially un- 



reduced ion-to-electron mass ratio of 16 . The effects of a guiding 
magnetic field have been cons idered by|N ishikawa et al. (2003 ); 
Hed edal & Nishikawa| ( |2005] ); |Dieckmann et al.| ( |2006| ) with 2D 
and 3D PIC simulations. 

GRB jets may carry with them a significant fraction of 
positrons (|Piran|1999]).|Kazimura et al.| (l998); J aroschek et al.| 
( 2004|); [Spirkovsky| ( |2008| ) modelled with 2D and 3D PIC sim- 
ulations the collision of two unmagnetized clouds, each con- 
sisting of electrons and positrons. Hoshin o~et al.| ([l992) intro- 
duced heavier ions in ID simulations. Magnetic field effects on 
such collisions were taken into account by |Spitkovsky| ( [2005] ); 
Sironi & S pitkovsky] ( |2009|). Interpenetratin g leptonic plasmas 



have been investigated by Silva et al. ( 2003 ). 

The simulation results typically show the generation of mag- 
netic fields by a filamentation of the plasma that implies a sep- 
aration of the currents, provided that the guiding magnetic field 
is not too strong and that the flow speeds are relativistic ( |Cary| 
et al.||198T] ). The energy density of the magnetic field reaches 
typically about 10% of the leptonic flow energy. However, most 
simulation studies could not observe a suprathermal population, 
as long as the plasma cloud collision speeds or the beam speeds 
are mildly relativistic. Such an acceleration of electrons on short 
spatiotemporal scales is conditional on both the presence of ions 
in the plasma flow and a mechanism that can transfer a signifi- 
cant fraction of the ion flow energy to the electrons. 

Initial conditions that result in a substantial magnetic field 
amplification and in the electron a cceleration to ultrarelativis - 
tic energies have been proposed by Bessho & Ohs awa| ( |1999| ), 
who studied in one spatial dimension the collision of magne- 
tized plasmas at the speed 0.9c and with a magnetic field direc- 
tion tilted by 45 degrees relative to the flow velocity vector. An 
ion-to-electron mass ratio of 100 was used. The acceleration of 
electrons up to Lorentz factors of -130 has been found. 

Here we consider the collision of two plasma clouds with 
a density ratio of 10, with a speed of 0.9c, and in the pres- 
ence of a strong magnetic field. The initial magnetic field would 



correspond to a primordial jet magnetic field (Lyuti kov et al. 
2003 ; Granot 200 3|. These initial conditions are similar to those 



magnetized plasma clouds colliding with a density ratio of 3, 
both composed of electrons and ions. The Lorentz factor of the 
collision speed was 2-3. Computational constraints demanded a 



used by |Bessho & Ohsawa ( [1999 ). Our magnetic field direction 
is, however, tilted at 0.1 radians relative to the flow direction 
and, most importantly, the two-dimensional simulation permits a 
more complex array of physical processes at a higher mass ratio. 
A ID study using almost the same parameters (Dieckm ann et al.| 
2008 ) suggests that the magnetic and electron energy density 
will increase drastically in the forming shock transition layer. 
This expectation is confirmed here, but we will demonstrate that 
the actual plasma dynamics in our 2D simulation differs notably 
from that in the previous ID studies. 

We demonstrate that the coherency of the circularly polar- 
ized electromagnetic wave ahead of the shock is reduced com- 
pared to what was assumed for the ID simulations, which can 
be partially attributed to the failure of the guiding magnetic field 
to suppress the filamentation. The shock planarity, which is en- 
forced by a ID simulation, is destroyed here by the development 
of a flux tube. This flux tube is the focus of the research let- 
ter [Murphy et al. ( 2010b ), hereafter MDD. Here we examine in 
more detail this flux tube and the plasma conditions in the PIC 
simulation that result in its growth. 

In this paper we perform a numerical simulation of two 
plasma clouds colliding, using conditions that are probably ap- 
propriate for GRBs. In Section [2] we describe the method used, 
in Section [3] we present some results and finally in Section]?] we 
discuss the results obtained. 
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Fig. 1. Sketch of simulation initial conditions 

2. The Numerical Experiments 

The particle-in-cell (PIC) simulation method has been described 
in detail elsewhere ( [Daw son 1983 ). The plasma is represented 
by an ensemble of computational particles (CPs). These CPs 
correspond to phase space blocks rather than physical particles. 
Consequently, the charge qt and the mass nti of a CP of the 
species i do not have to correspond to the equivalent of the phys- 
ical particles it represents. However, the charge to mass ratio 
qi/rrii must equal that of the physical particle. Then the ensemble 
properties of all CPs representing the species i are an approxi- 
mation to the ensemble properties of the corresponding plasma 
species, e.g. an electron species or an ion species. 

PIC code s approximate the Klimontovich-Dupree equations 
( Dupre ep963| ), which correspond to the solution of the Vlasov- 
Maxwell equations by the method of characteristics. PIC codes 
can capture all kinetic wave modes and instabilities found in col- 
lisionless plasma, if the simulation box size and the resolution, 
as well as the statistical plasma representation, are adequate. 

We consider the collision of a dense plasma cloud with a ten- 
uous one. The electrons with mass m e and the ions with mass 
nti = 250m e of the dense cloud are the species 1 and 2, re- 
spectively. The electrons of the tenuous cloud are species 3 and 
species 4 denotes the ions of the tenuous cloud. We normalise 
our variables with the plasma frequency of the species 2, with 
the density ri2, the charge and the mass mi. The normaliza- 
tion is useful, in that it renders the simulation results indepen- 
dent of the plasma density, which is unknown for GRB jets. The 
skin depth of species 2 is \i = c/a) P 2. The elementary charge 
is e. The quantities in SI units (subscript p) can be obtained by 
the substitutions E p = oj^cm^je, = co^mfil e, p p - en^p, 
J p = ecn^S, x p = A2X and t p = t/a) P 2. The solved equations are 



V x E = -d t B, V x B = d t E + J, 
V • B = 0, V • E = p, 



dtVj = qj (E + Yj X B) , Pj - = mjTvj, d t Xj = v 7 - 



(1) 

(2) 

(3) 



Here the subscript j refers to the f h CP with the mass mj 
and the charge qj. We use the Plasma Simulation Code PSC for 
our simulations, which is a relativistic 3d MPI parallel domain 
decomposed PIC code. It has been extensively used in the laser 
plasma community ( [Roth et al.|2001[ [Cowan et aL|2 004). 



2.1. Simulation Initial Conditions 

We begin the simulation with two colliding plasma clouds. 
Figure [T] illustrates the flow and field geometry. The inclina- 



tion angle relative to the flow velocity vector of the initial mag- 
netic field Bo is 0.1 radians and its magnitude |Bo| = R l/2 in 
our normalization, where R = mi/m e = 250 is the reduced 
ion/electron mass ratio. The magnitude |Bo| is thus such, that 
the electron cyclotron frequency equals the electron plasma fre- 
quency of the dense plasma cloud. Both beams (high and low 
density) travel initially at 0.63c but in opposite directions, giv- 
ing a relative speed of 0.9c. This speed jump will be distributed 
over the forward and reverse shocks. The temperature of all 
species is T = 131keV. The thermal velocity of the electrons is 
v t h,e = ^JkT/m e = 0.83v£ = 0.52c. The thermal velocity is such 
that it is not greatly smaller than the collision speed, as expected 
to be for internal shocks in GRB jets. The thermal velocity of 
the ions is v t h,i = Vth,el V# = 0.033c. The initial distribution of 
the particles is a relativistic Maxwellian or Maxwell- Juttner dis- 
tribution. The initial density ratio is chosen to be 10. We cannot 
use here the piston method of |Forslund & Shonk| ( |1970| ). In this 
piston method, the plasma is reflected at a conducting wall. The 
plasma symmetry across the wall is exploited and only one cloud 
has to be modelled. This cloud is reflected onto itself by the wall 
and a shock develops at the collision boundary. It is a compu- 
tationally efficient method. However, only a limited number of 
field geometries should be modeled with it. The electric fields at 
a conducting wall must point along the surface normal and the 
magnetic field orthogonal to it. We can also model with it un- 
magnetized plasma as in |Forslund & Shonk (1970]). It is further- 
more not possible to let an oblique magnetic field stream with 
the plasma into the simulation box through an open boundary 
and let the electromagnetic field distribution at the conducting 
wall develop self-consistently when the plasma reaches it. The 
oblique geometry implies a flow-aligned magnetic field com- 
ponent, which would have to end at the front of the inflowing 
plasma cloud. A magnetic monopole would be the consequence. 
While the particular choice of the boundary condition may not 
affect the long term evolution of the shock far from the boundary, 
a representation of both colliding plasma clouds will be more 
realistic. In our case study, we collide two plasma clouds with 
different densities and, in fact, this asymmetry rules out the pis- 
ton method altogether due to an absent mirror symmetry at the 
wall. This is physically motivated by the expectation that plasma 
clouds of similar but unequal density will collide in the GRB 
jet, although we have to point out that the collision boundary is 
unlikely to be as abrupt as the computationally convenient one, 
which we implement here. However, we do not expect a strong 
dependence of the key simulatio n results on the boun dary shape. 
The smooth boundary used by |Bessho & Ohsawa (1999) and 
the sharp one used by |Dieckmann et al.| ( |2008| ) did not result 
in qualitatively different simulation results. The simulation will 
furthermore show that the important structures form well after 
the initial time, when the boundary has been smeared out. The 
transport of the magnetic field B at velocity u || x gives a convec- 
tion electric field E c = u x B (Baumjohann & Treu mann|1996| ). 
It changes its direction at the cloud collision boundary. 



2.2. Linear instability 

The original motivation to pick the high initial plasma temper- 
ature, the strong guiding magnetic field B and the asymmetric 
beam density has been to suppress the filamentation instab ility 
and, thus, to enforce a planar shock (Dieckmann et al. 2008 ). We 
can test this assumption by computing the approximate linear 
dispersion relation under the following assumptions. An early 
time is considered, when both clouds overlap in a small interval 
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is resolved in 1 cell in the simulation. 




Fig. 2. The exponential growth rate of the linear instability in 
units of co P 2 for the ion-to-electron mass ratio R = 250 used 
in the simulation. The wavenumbers are normalized to the skin 
depth A2 of the dense ion species. 



along x. The interval is large enough to ensure that the elec- 
trons of both clouds have mixed and form a single, spatially 
uniform and hot distribution. The interval is small enough so 
that firstly the ion distribution is unchanged and secondly the 
magnetic field component orthogonal to the beam velocity vec- 
tor has not yet been compressed to a significant amplitude. Then 
we can approximate the plasma in the cloud overlap layer by 
two counter-propagating ion beams, which move through a hot 
electron background along a guiding magnetic field with an am- 
plitude that equals 250 1/2 . The relative speed between both ion 
beams is 0.9c and their density ratio is 10. 

Figure|2]displays the linear dispersion relation calculated for 
the ion-to-electron mass ratio of R = 250 in the simulation. 
The exponential growth rates peak in the field-aligned direction, 
which is characterized by a wavenumber component along the 
beam velocity vector k\\A.2 ± and a perpendicular component 
k ± A2 = 0. These modes could be observed in a 2D simulation 
of an oblique shock, which employed a lower collision speed 
( [Dieckmann et al.|2010| ). However, the modes with & ± /i 2 ^ are 
not suppressed. Bands of unstable waves reach out to k ± A,2 » 1. 
The shock and its downstream region may thus not be planar. 
According to this solution of an idealized linear dispersion rela- 
tion, the structuring of the shock along its boundary is captured 
well by a simulation box that spans a few ion skin depths into 
this direction and resolves the electron skin depth. Our simula- 
tion will represent the wavenumber band 1 ~ k ± A>2 < 134. The 
growth rate map for the mass ratio R is qualitatively similar to 
that obtained for the correct proton-to-electron mass ratio (not 
shown), which suggests that the spectrum of unstable waves in 
the PIC simulation will be realistic, at least during the initial 
time. This is not always the case for reduced mass ratios ( [Bret &| 
|Dieckmann|20T0l ). 

2.3. Numerical resolution and computational details 

For the 2D simulation, the box measured in ion skin depths 
is of width L x = 656/12 and of height L y = 6A2 resolved in 
2.8 x 10 4 cells in the propagation direction and 256 cells in 
the perpendicular direction. The total plasma Debye length L T , 



where j- = jj— + 77- 

T L D,ion L D,electron 

The electron skin depth is resolved in 2.7 cells. We use 200 com- 
putational particles per cell (100 ions, 100 electrons) in the dense 
plasma and 100 particle per cell (50 ions, 50 electrons) in the ten- 
uous plasma, which is possible by assigning different numerical 
weights to the CPs. No new particles are introduced at the pe- 
riodic boundaries during the simulation. The two plasma clouds 
rapidly detach from the boundaries. The high thermal speed of 
the electrons implies, that they leave the plasma clouds at their 
rear ends to leave behind a positive net charge. This charge sep- 
aration induces an electric field, which accelerates the ions. This 
process has been researched in the context of a plasma expansion 
into a vacuum ( Mora & Grismayer 2009 ). While such an expan- 
sion clearly is an artifact of our initial conditions, it does not 
visibly influence the plasma dynamics at the collision boundary. 
This will be demonstrated below by the supplementary movies, 
which do not show waves or plasma structures that propagate 
from the simulation box boundaries to the cloud collision bound- 
ary in the center of the simulation box. The x-boundaries are in- 
tentionally placed sufficiently far from the shock forming region 
that no signal can reach the shock forming region traveling from 
the boundaries within the simulation time. 

The number of processors used was 256 Intel Xeon E5462 
2.8GHz on an ICHEC SGI ICE machine (Stokes). Total wall 
clock time was 156 hours, giving a cumulative CPU time of ~ 
40, 000 hours for the 2D simulation. 



3. The simulation results 

In this section we elucidate the consequences of two plasma 
clouds colliding. Although the plasma dynamics in the cloud 
overlap layer is determined by both clouds simultaneously, we 
will analyse them separately. This is made possible by tagging 
the CPs initially belonging to the species under consideration. 
Henceforth we shall refer with dense ions/electrons to the cor- 
responding species of the dense cloud that moves to the right, 
while the tenuous ions/electrons are those of the diluted plasma 
cloud that moves to the left. In order to gain insight into the 
growth of structures due to the interaction of the clouds, we 
present results of a numerical simulation at the early time of 
Tl = 62 and the later time of T2 = 182. The static images are 
supported by animated MPEGs available in the online edition, 
which show the time-evolution of the fields. 



3.1. Early stage at time T1 

At the time Tl, the beams have counterstreamed for tens of ion 
skin depths. We look first at the electron phase space distribution 
integrated over y. The dynamics of the electron beams show that 
the electrons at position x=30 are already accelerated to T = 60 
from the initial T < 8 (Figure [3]). The strongest electron accel- 
eration occurs for 30 < x < 33. The relativistic inertia of these 
electrons is already comparable to that of the ions and we expect 
a detectable reaction from them. From the ion phase space plot 
in Fig. |4| which is also integrated over y, it can be seen that lit- 
tle interaction has taken place in between the colliding clouds of 
ions over most of the displayed x-interval. However modulations 
of the specific x-momentum of both clouds are highly evident at 
30 < x < 40. All ions are decelerated at 30 < x < 33 in the 
simulation frame, which provides a reservoir of energy for the 
acceleration of electrons and dense ions to Tv x « 1.2 at x « 32. 
The dense ions are decelerated also at x > 35. 
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3. Electron phase space: Log of electron density as a func- 
of Lorentz factor T and x at time t-T\. 




Fig. 4. Ion phase space: Log of ion density as a function of spe- 
cific x-momentum Tv x and x at time t = T\ 



The large-scale distribution of the ion densities is shown in 
Fig. [5] for the tenuous ions and dense ions. The tenuous ions 
undergo a rapid filamentation, as they reach the cloud overlap 
layer at x « 40. The tenuous ions initially form filaments with 
a thickness of « 0.1 that are almost aligned with the flow veloc- 
ity vector. Assuming the current channels are engendered by a 
filamentation instability between the ions of both clouds with- 
out electron involvement, one should see such structures also in 
the dense ions. This is because both clouds have approximately 
the same density and temperature for 30 < x < 40 and both 
should thus behave similarly. Filamentary structures resembling 
those in the tenuous ions are, however, not visible in the dense 
ions on this scale. This filamentation instability must thus in- 
volve all plasma species. The electrons in this interval still carry 
a substantial directed flow energy, which can be released and fur- 
ther modify the instability. At this advanced simulation time, the 
linear dispersion relation discussed above would no longer be a 
good approximation. However, the key conclusion we draw from 
it, namely that the filamentation instability is not suppressed, is 
supported by the ion distribution. 

The current channels in the tenuous ions are then rapidly de- 
flected and thermalized, as they enter the interval 30 < x < 33 
in Fig.[5j where the electron acceleration and, thus, the electro- 
magnetic fields are strongest. A structured and beamed distri- 
bution is present for x < 30. The ions in these channels have 
not decreased significantly their p x momentum (Fig. [4]). A two- 
dimensional density modulation is also visible in the dense ions 
for -20 < x < 30. Its density peaks at x = 30, which is the 
interval where the dense ions are slowed down most in Fig. [4] 
The planar front of the dense cloud is surrounded by filamentary 
distrib utions, as it has also bee n observed for a reduced collision 
speed (Dieckmann et al.|2010| . 




Fig. 5. Two dimensional logarithm of the ion density as a func- 
tion of X and Y at the time t = T\ . The upper panel shows the left 
moving tenuous ions and the lower one the right-moving dense 
ones. 



A zoom of the zone around the collision boundary (Figure [6]) 
shows the spatial density distribution of each of the four species. 
The current channels in the tenuous ions start to form at x « 38. 
The current density increases as we go from this position to x « 

34 and the space between the channels is progressively depleted 
of ions. The current channels merge, e.g. at x « 33 and y « 
4. Only 4 major current channels eventually cross the electron 
acceleration region and reach x « 30, where they are deflected 
to increasing values of y and scattered at x « 27. Remarkable 
density modulations of the tenuous ions are visible at x « 30, 
where we find a minimum density of exp(1.3) « 4 at y « 1.7 
and a maximum density of exp(3.9) ~ 50 at y « 2.5. Their 
density minima at (x,y) = (29.5, 1.7), (30,3) and (30,5.5) are 
correlated with local maxima of the density of the dense ions. 
Some correlations between the filaments of the tenuous and the 
dense ions are visible, e.g. at x « 37 and y « 1 and y « 2, where 
the current channels both have a density of « exp (3.5). 

The dense ions show a quasi-planar front at 30 < x < 3 1 with 
an enhanced density as well as a less pronounced second front at 
x « 36. In particular the front at x « 31 is sharp and the densities 
of the dense ions and electrons decrease by a factor « 5 and « 20, 
respectively, as we cross its boundary to the upstream (larger x). 
The lower density of the dense electrons ahead of the front is 
probably caused by their rapid expansion upstream. Figure [^re- 
veals that hot electrons leak out of the acceleration region and 
reach a x « 50. Their current induces a return current in the ten- 
uous electrons, which accelerates the latter towards the density 
peak at x = 31 and results in their accumulation at 30 < x < 33 
and at 35 < x < 38. The dense electrons also accumulate in 

35 < x < 38. The density maxima of the tenuous electrons cor- 
relate well with the locations, where the electron acceleration in 
Fig. [3] and the ion deceleration in Fig. [4] is strongest. The ion 
deceleration in the interval 35 < x < 38 causes their accumu- 
lation and also that of the electrons, which must cancel the ion 
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Fig. 6. Two dimensional plots of logarithm of dense and tenu- 
ous ion and electron clouds at time Tl, for the zone close to the 
collision boundary. 



charge. Furthermore, a strong filamentation along y is visible for 
20 < x < 30 for all four species in Fig. [6] 

The mechanisms that accelerate the electrons at the expense 
of the ion energy can be identified with the help of the electro- 
magnetic fields. The magnetic and electric components in the 
location of the collision boundary (Figures. [7] and [8]) reveal a 
large electromagnetic pulse in the interval 30 < x < 33, which 
is the interval with the strongest electron acceleration and ion 
deceleration in Figs. [3] and [4] The electromagnetic fields orthog- 
onal to the collision (x) direction reveal bipolar pulses, which are 
shifted in space. This is most evident for B y , where the negative 
and positive magnetic field is separated at x « 31.5. This zero 



crossing of B y is approximately where B z reaches its maximum 
positive value, evidencing a phase shift of 90° along x between 
By and B z . A negative B z is then visible at x « 33, just ahead of 
the positive interval of B y . The B y and B z reach moduli ^100 
each, which is about 6 times larger than |Bo| = 250 1/2 . 

The bipolar pulse of E y is in phase with that of B z , while 
the pulse in E z is in antiphase with that in B y . The electric field 
pulse is to some degree the consequence of the rapidly moving 
magnetic field pulse (convection electric field) and this contribu- 
tion would vanish in the rest frame of the magnetic field struc- 
ture moving with the speed « v h . The B y and E z show a further 
oscillation at x « 35, suggesting that it is the same circularly 
polarized energetic electromagnetic structure that was observed 
previously byjDieckman n et al. ( 2008| 2010| ), although here its 
coherency along y is low and it does not extend far upstream. The 
low coherency can at least in parts be attributed to the strong spa- 
tially varying B x . The presence of a strong nonplanar B x implies 
spatially varying strong currents in the z-direction. 

A single circularly polarized purely electromagnetic wave 
cannot by itself accelerate electrons to ultrarelativistic speeds. 
Its combination with the electric E x » in Fig. [8 is necessary 
for this purpose. In our normalization the ratio of the values for 
E and B is almost equal to the force ratio for a particle, which 
moves with v « c. The electric and magnetic forces on relativis- 
tic particles are thus comparable. The E x > can be explained 
as follows. As the upstream plasma impacts on the strong mag- 
netic field at x « 33, the electrons are deflected away from their 
original flow direction, while the ion reaction is much weaker. 
A current develops in the y-z plane, which amplifies locally the 
magnetic field. The electrons fall behind the ions, because their 
velocity along x is reduced. An E x > builds up, which tries to 
restore the quasi-neutrality. The tenuous electrons are dragged 
by it across the magnetic field. This cross-field transport accel- 
erates the electrons to relativistic speeds, as it is confirmed in 
Fig. [3] This acceleration acts in the y-z plane, which partially 
explains the spatial confinement along x of the accelerated elec- 
trons. This spatial confinement of the strong electromagnetic 
fields and of the electron acceleration also implies a localized 
fast decrease of the ion flow speed, explaining the ion reaction 
within 31 < x < 35 in Fig. [4] The electric field, which drags 
the electrons to the left in this interval, causes the slowdown of 
the tenuous ions as they move to the left and the speedup of the 
dense ions, which move to the right. The consequent ion accu- 
mulation further confines with its massive positive charge the 
electrons. This electron acceleration mechanism ceases to work, 
when the E x is strong enough to stop the ions. 

We can also observe filaments in particular in B z for x < 30, 
which results out of the filamentary strucures observed in all 
plasma species in this interval. We expect this magnetic compo- 
nent to show the strongest modulation if the currents are approx- 
imately aligned with the x-direction and are modulated along 
y. Magnetic field stripes aligned with x are also observed for 
x > 35. These filaments are driven by the upstream electrons 
and the electrons that leak to higher x, reaching x « 50 in Fig. [3] 
This upstream filamentatio n has a lso been observed at the faster 
shock modelled by jMartins et al.| ( [2009] ). 



3.2. Late stage 

The physics at time Tl has given us insight into the plasma pro- 
cesses that take place initially and precondition the plasma such 
that the vortex forms, which is discussed in detail in MDD| We 
now consider the simulation time T2, which is approximately 
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Fig. 7. Magnetic field components: Zoom of 2D linear Fig. 8. Electric field components: Zoom of 2D linear colourscale 
colourscale plot of spatial distribution of B x (upper plot), B y plot of spatial distribution of E x (upper plot), E y (middle plot) 



(middle plot) and B z (lower plot) at t = T\. 



and E z (lower plot) at t = T\ . 



Ion Phase Space t= 182.4 



that investigated in MDD Here we discuss in more detail the 
electromagnetic fields at the shock front, which have not been 
shown in MDD The fields may not provide more information 
related to the vortex than the current, but they are essential to 
understand the particle acceleration, which is the focus of this 
paper. 

To aid the analysis, in the simulation three distinct zones are 
identifiable at T2. These zones are most easily distinguished in 
Fig. [9] which displays (as in Fig. [4} the ion phase space den- 
sity integrated over y as a function of x and p x , but now at the 
time T 2 . The displayed x-interval can be subdivided into the in- 
terpenetrating ion beam zone (IIBZ) with -35 < x < 60, the 
downstream region with 60 < x < 80 and the foreshock region 
of the strong forward shock with 80 < x < 130. Fig. [4] and Fig. 
TT are show a similar time slice to Fig. 2 in MDD| 

First we note that the strongest interaction takes place at 
x » and that it is tied to the front end of the dense cloud. 
No energetic structure is visible at the front end of the tenu- 
ous cloud (not shown). The incoming plasma from the upstream, 
here the tenuous cloud, is reflected at x « 80, forming a shock- 
reflected ion beam. This shock reflected ion beam is hot and the 
reflection is not specular. A bunch of ions with values of p x , 
which are comparable to the initial ones of the dense ions, is vis- 

















50 100 



Fig. 9. Ion phase space: Logarithm of electron density as a func- 
tion of specific x-momentum Tv x and x at t = 7^. 



ible in the interval 100 < x < 120. These are the dense ions, 
which were located ahead of the strong interaction region found 
at x « 33 in Fig. [4] and they have not been accelerated by it. 
The downstream region is characterized by a single, almost spa- 
tially uniform and hot ion population. The kinetic energy stored 
in the relative speed between the upstream and the downstream 
plasma has been converted into heat. The momentum conserva- 
tion, together with the asymmetric cloud densities, implies that 
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Fig. 10. Ion momentum distribution at x = 73.2, normalized to 
its peak value, at t = T 2 . A fit to a Maxwell- Jiittner distribu- 
tion with the thermal speed 0.18c and the mean speed 0.33c is 
overplotted. 



the downstream region cannot be stationary in the box frame. 
Indeed, the normalized mean speed along x of the downstream 
ions is « 0.33 according to Fig. [10| while the speed modulus of 
both incoming clouds is the same in the simulation frame of ref- 
erence. The speeds of the observed forward shock that is moving 
to the right and the reverse shock that is still developing, which 
are given by the relative speed between the downstream plasma 
and the respective upstream plasma, must differ. 

The IIBZ to the left is characterized by the co-existence of 
the dense ions and the tenuous ions, which have crossed the front 
of the dense cloud prior to the formation of the shock. The mean 
speed modulus of the tenuous ions increases as we go to lower x 
and it is close to the initial one at x « -35. This spatially vary- 
ing mean momentum is a consequence of the shock formation. 
The mean speed modulus of the tenuous ions in the strong inter- 
action region decreased steadily in time, as the electron acceler- 
ation became more efficient. The later in time the tenuous ions 
traversed this strong interaction region, the more energy they lost 
to the electrons. Once the downstream region has been formed, 
the tenuous ions can no longer cross this obstacle. The tenuous 
ion beam in the IIBZ is thus a transient effect. Eventually, a re- 
verse shock will form between the IIBZ and the downstream, 
giving rise to a shock-reflected ion beam. 

Movie 1 shows a zoom of the time evolution of the ion dis- 
tribution in (p x , x) space. We note that no signal, either wave 
or plasma structure is detected propagating into the box either 
from left or right boundaries. The early stage shows two beams 
colliding. The beams interpenetrate and then decelerate and a 
forward shock begins to form. The tenuous beam is partially re- 
flected to high velocities by the dense beam. Left of the collision 
boundary, the structure is clearly heated and takes on a thermal 
distribution (fitted to a Maxwellian in Fig. [lOb. The lower panel 
shows the left moving cloud ion distribution in (x, y) space. The 
most significant structures at early times are the distinctive fila- 
ments, which dominate both the foreshock and the downstream 
region. The filaments are sheared in the negative y direction and 
circular structures are seen to form which increase in size at the 
later stage of the simulation, to eventually fill the simulation box, 
at which point we stop the simulation. 

Figure [TT] displays the electron phase space density distribu- 
tion as a function of x and T at T2. The downstream region hosts 




Fig. 11. Electron phase space distribution at t = T 2 : Logarithm 
of electron density as a function of the Lorentz factor T and x. 



a hot and dense electron population. The strongest electron ac- 
celeration takes place at the boundary between the downstream 
region and the foreshock at x « 80. The most energetic elec- 
trons reach T > 180. Their relativistic energies are well above 
that of an ion with 250 electron masses and the speed v&. The 
latter would have a kinetic energy, which would equal that of an 
electron with T « 70. A dilute cloud of ultrarelativistic electrons 
is leaking out from the downstream region into the foreshock re- 
gion and into the IIBZ region. The distribution of the dense elec- 
trons in the IIBZ is practically unchanged and clearly separated 
from the hot population, indicating that they are not interacting 
strongly through beam instabilities. The turbulent tenuous elec- 
trons in the foreshock with x > 80 have probably been heated by 
the shock-reflected ions, since they are perturbed in the spatial 
interval up to x « 130, which coincides with the cut-off of the 
ion beam. Some interaction may have taken place in form of a 
filamentation instability between the leaking hot electrons and 
the incoming tenuous electrons, as it was observed at the time 
Tl. 

Movie 2 shows a zoom of the time evolution of the electron 
(r, x) phase space in the upper panel, and the distribution of the 
left moving cloud in the lower panel. We note that no plasma 
structure enters the box from the left or right x- boundaries. The 
acceleration to ultrarelativistic Lorentz factors occurs in a con- 
fined region close to the collision boundary. Filamentation dom- 
inates but to a lesser extent than found in the ion clouds, due to 
the higher thermal velocities of the electrons in the foreshock. 
At later times, the same circular structure is seen as in the ion 
clouds. 

The magnetic energy density at the time T2 is displayed in 
Fig. 12 Elevated magnetic energy densities are observed in the 
IIBZ up to x « 50. Then the magnetic fields strengthen and reach 



a magnetic energy density plateau with 20|B | at x « 60, which 
spans the downstream region up to x « 70. Then a massive 
peak is observed, reaching a peak value of 10 2 |B | 2 at x « 73. 
The magnetic energy density decreases rapidly as we go from 
this position to increasing x and it reaches a local minimum at 
x « 85 with the value « 2|B | 2 . The second weak maximum at 
x « 90 is followed by an apparently exponentially decreasing 
magnetic energy density that cannot be distinguished from |B | 2 
at x « 125. The absolute maximum of the magnetic energy den- 
sity is found a few ion skin depths to the left of the x « 77, in 
which the upstream ions are reflected in Fig. [9] and where the 
electrons experience their strongest acceleration in Fig. [TT] It is 
thus the magnetic ramp in Fig. 12 that is responsible for the par- 
ticle acceleration, but the plasma structure responsible for the 
extreme magnetic energy density is located behind it. The po- 
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Fig. 12. The total magnetic energy density B 2 at time T2 aver- 
aged over y and normalised with respect to |B | 2 plotted against 
x. The field is amplified to above 50 times its initial value in a 
thin layer, far above the value expected from a shock compres- 
sion. 
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Fig. 13. Two dimensional logarithm of the ion densities. The 
upper panel shows the left-moving tenuous ions and the lower 
panel the right-moving dense ions at t = 7^. 



sition x sb 125, beyond which the magnetic field is not visi- 
bly amplified, is approximately colocated with the front of the 
shock-reflected ion beam in Fig. [9] 

sug- 



The large magnetic energy density observed in Fig. 12 



gests, that the underlying currents must be due to the ions and 
we analyse their density distribution in Fig. 13 A structure is ob- 
served at-10<x<0 and y < 2, which is periodically wrapped 
around at the boundary y = 0. At the same position, we find a 
weak magnetic energy density peak in Fig. [12] A strong density 
depletion is only seen in the dense ions. The movie shows that 
this structure can be interpreted as the two-dimensional equiva- 
lent of a magnetic bubble. The interplay of the current filaments 
results in the accumulation of magnetic energy in a localized 
pocket, which is convected with the dense ions. The pressure 
gradient force of this magnetic bubble expels the dense ions. The 
tenuous ions move with a relativistic speed in the rest frame of 



this bubble and they and their density distribution are practically 
unaffected by the magnetic pressure gradient force. 

The interval 50 < x < 73 with the extreme magnetic en- 
ergy is characterized by the increased ion density, which we ex- 
pect to find in the downstream region of shocks. Note that the 
shock compression is not unusually high. The ratio between the 
downstream ion density sb exp (5) and the upstream ion density 
sb exp (4) (summed over both ion species) at the forming reverse 
shock at x sb 40 is 3 in Fig. [13] The ratio between the down- 



stream density and the upstream ion density at the forward shock 
at x w 75 is about 4. No current filaments are observed in this in- 
terval and it is unclear from this plot what is behind this immense 
magnetic energy density. 

The tenuous ions show filamentary structures in the interval 
80 < x < 130, in which we find the weaker peak in the mag- 
netic energy density in Fig. [12] This peak is thus due to the ion 
beam filamentation instability. The narrow well- separated fila- 
ments, which are stretched out over tens of ion skin depths, can 
be explained with a quasi-equil ibrium similar to that derive d for 
relativistic electron beams by Hammer & Rostocker ( 1970). No 
filaments can be found in the dense ions in this interval. The 
dense ions contribute to the fast beam, which outruns the for- 
ward shock, and this beam is too hot to react to the magnetic 
fields. The seemingly exponential growth in space of the mag- 
netic energy density in Fig. 12 as we go from x ^ 130 to x sb 95 
may result from the combination of the exponential growth in 
time of the ion current in response of the filamentation instabil- 
ity with the convection to the left of the tenuous cloud. 

Now we turn to a more thorough examination of the struc- 
ture, which is responsible for the extremely strong magnetic 
field. It has previously been identified in MDD as a flux tube. 
Figure 14 displays the three magnetic field components at the 
time T2. The B x and the B y magnetic fields form a magnetic loop 
in the x-y plane for 70 < x < 77, while the B z component reveals 
an elliptical structure within 71 < x < 76. The simulation geom- 
etry implies that we can understand this magnetic field geometry 
as the combination of the axial magnetic field of a coil with an 
infinite extent along z and a magnetic ring that is surrounding 
it. Such flux tube distributions can be force-free, as they fullfill 
V x B = aB with a constant a. If such a structure were to form 
in an astrophysical jet far from any magnetic source object, like 
what remains from the progenitor star, the magnetic field lines 
in the z-direction would have to be closed. A simple geometry 
that fullfills this necessary closure is a spheromak that resembles 
a smoke ring. 

Apart from the dominant flux tube distribution, the mag- 
netic field shows further elongated structures in the interval 
55 < x < 70. The geometry of these structures is reminiscent of 
the end product of the filamentation instability, which is appar- 
ently thermalizing the downstream distribution while maintain- 
ing the strong magnetic fields in this domain in Fig. [12] However, 
once the plasma has thermalized, these magnetic fields can no 
longer be upheld by plasma currents (Waxman 2006). 

The strongest electron acceleration and the ion reflection by 
the forward shock takes place at x sb 77, where we find mag- 
netic stripes in the B y and B z components. These structures are 
not closed and resemble the bipolar pulse we have observed at 
the earlier time Tl. If this bipolar pulse is responsible for the 
particle acceleration also at this late time, we expect to find a 
strong positive electric E x component at x sb 77. Figure [15] < 



con- 



firms this. While the bipolar pulse is responsible for the electron 
acceleration and for the current sheet, out of which the flux tube 
has grown, the flux tube itself scatters through its magnetic field 
the electrons, as we can see from Fig. 1 1 



10 



G. C. Murphy et al.: Field amplification and electron acceleration in a protoshock 



B If t= 182co ni 



5 E- 

4 i- 



v» * 




60 65 70 75 



60 65 70 75 




82.6 128.0 



-108.6 -51.4 5.8 62.9 120.1 177.3 




32.5 -16.5 -0.4 15.7 31.8 47.8 



Fig. 14. Magnetic field components at the time T2: Zoom of a Fig. 15. Electric field components at the time T2: Zoom of a 
2D linear colorscale plot of the spatial distribution of B x (upper 2D linear colorscale plot of the spatial distribution of E x (upper 



panel), B y (middle panel) and B z (lower panel). 



panel), E y (middle panel) and E z (lower panel). 



The separation of the flux tube from the current sheet driv- 
ing it may have an important consequence. Figure [14] demon- 
strates that the periodic boundary conditions limit the growth of 
the flux tube. Selecting a simulation box that is much larger in 
the v-direction may result in a larger flux tube, but only if the fur- 
ther growth is not limited by the thickness of the current sheet 
along x. We would expect such a size limitation, if the flux tube 
could only exist in the current sheet. The current sheet would, of 
course, widen, if the simulation would model protons rather than 
the lighter ions. The thickness of the current sheet would, how- 
ever, still be limited by the distance, over which the ion energy is 
depleted by the accelerating electrons. Here the simulation sug- 
gests that a current sheet ahead of the flux tube suffices to drive 
it. The flux tube can thus grow to a large MHD size and be a 
reservoir of magnetic energy that may not be dissipated away as 
quickly as that due to current filaments. 

The electric field furthermore demonstrates strong electric 
fields, which are partially correlated with the magnetic fields of 
the large flux tube and of the downstream filaments. The E x and 
E y components show, for example, the same topology as the flux 
tube's B z component, while the E z component resembles the flux 
tube's By distribution. The electric field amplitude is well below 
that of the magnetic field. Those of E x and E y are about half 



that of their magnetic counterparts, while that of E z amounts to 
a third of B z . The electric energy density will thus be about 20% 
of the magnetic one. 

Movie 3 shows a zoom of the time evolution of all six 
components of the electric and magnetic fields. We note that 
no signal, either wave or plasma structure, can be seen enter- 
ing the simulation box on the left or right x-boundaries. In the 
lower four panels, E y , E z are clearly seen to be in phase and 
antiphase, respectively with the B z ,B y components. Initially a 
strong dipolar field dominates, and filaments are clearly visi- 
ble in the E y ,B z components. The components parallel to the 
shock normal (B x , E x ) are considerably more fragmented than 
their perpendicular counterparts. 

3.3. Generation of Mid-Infrared Synchrotron emission 

The relativistic electrons are expected to emit synchrotron emis- 
sion in the presence of a magnetic field. The synchrotron roll- 
off frequency is at v ro iioff = ^c7 2 > where oj c is the electron cy- 
clotron frequency. In the simulation we have found that for a 
mass ratio of R = 250 the electrons are accelerated to Lorentz 
factors of ~ 200. We may assume from a comparison of the 
ID PIC simulations with a mass ratio 100 and 400 in Refs. 
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IBessho & Qhsawal ([1999]) and |Dieckmann etaLl ( [2008] ) that the 
peak electron Lorentz factor is proportional to the mass ratio. In 
this case and for a full mass ratio many electrons would reach 
Lorentz factors of 1500, this gives 1500 2 6l> c . The magnetic field 
strength within a GRB jet is unknown and we have to resort to 
a guess. For an oj c = (jl> p , electron ~ 10 5 Hz this gives « 200 GHz 
in the rest frame of the jet. There is a relativistic Doppler shift 
(8 = V(l + P)I(\ -P)) from the moving frame in the jet into the 
observer's (Earth) frame. Assuming a/? = 0.99999 for GRB jets, 
giving 8 = 450, we derive a frequency of 100 THz in the infrared 
frequency range. Secondary processes such as inverse Compton 
scattering from the population of T ~ 400 electrons are expected 
to upscatter the photons to gamma ray energies. 

4. Discussion 

In this paper we have examined the collision of two plasma 
clouds at a shock speed of 0.9c. The aim of the study is to gain in- 
sight into the behaviour of mildly relativistic shocks, and the as- 
sociated phenomena of electron acceleration, magnetic field am- 
plification, filament formation. The shock speed is at the lower 
end of the interval proposed for internal GRB shocks (Piran 
|1999| ). We add a quasi-parallel guiding magnetic field to the sim- 
ulation, in order to probe the effects of a strong magnetization 
expected to be present. The mass ratio has been reduced to 250 
to make the simulation more computationally tractable, while re- 
taining the physical mass asymmetry which ensures a reservoir 
of ion energy is availa ble to contribute to particle acceleration 
( |Amato &~A rons 2006). The density ratio of 10 has been chosen 
to emphasize the effects of asymmetry in GRB shocks. This pa- 
per extends MDD discussing in more detail the conditions that 
result in vortex formation also examining the electromagnetic 
field distribution at the shock front. While MDD only considered 
the currents, in this paper we focus on the particle acceleration 
mechanisms, not covered in the previous publication. 

The magnetized background together with the colliding 
clouds ensures a jump in the convection electric field will be 
present. The initial jump in the convection electric field triggers 
the growth of a ramp in the magnetic field from which particles 
can be accelerated. This electromagnetic reaction to the plasma 
cloud collision is the first and most rapid response by the sys- 
tem. This ramp is initially planar as found by Dieckmann et al. 
( 2010| ), but later evolves through plasma instabilities into a non- 
planar structure. Previous studies have shown that magnetized 
shocks have little particle acc eleration unless ions are present 
( [Hoshino et al.|1992||Amato & Arons|2006| ). Equally studies of 
weakly magnetized collisions have shown little or no difference 
to unmagnetized plasma (Nishik awa et al.|[2QQ3j ). The simulta- 
neous presence of ions and of a strong guiding magnetic field 
show the rapid formation of a shock, which accelerates electrons 
to highly relativistic speeds and amplifies substantially the initial 
magnetic field. 

The results shown here - in the context of the observations 
of po larisation in GRBs ( [Coburn & Bog gs 2003 ; Steele et aL| 
2009) provide a compelling argument for the role of a dynam- 



ically significant primordial magnetic field from the jet. Such 
fields, in oblique collisions, provide a mechanism to transfer en- 
ergy down the mass scales from ions to electrons, allowing elec- 
trons to increase their relativistic mass until they can be injected 
into the Fermi acceleration mechanism. 

Based on the linear dispersion relation, we hypothesise that 
the filamentation instability is not suppressed in this regime of 
parameter space. The nonlinear simulation supports this hypoth- 
esis, but goes further to demonstrate a possible solution to the 



longstanding problem of filament lifetimes. Evidence for vortex 
formation is found, possibly a stable solution to offset rapid fila- 
ment decay mechanisms ( Waxman 2006). Such vortices are ob- 
served in nonrelativistic plasma flows ( Alexandrova et al.|2006 ) 
and may thus not be unlikely structures in the more energetic as- 
trophysical flows. The oblique shock considered here, also not 
unlikely in the context of GRB jets with helical background 
fields, allows sufficient transverse motion and transverse currents 
for a stable vortex structure to form through secondary instabili- 
ties. The simulation study reveals that the vortex has an internal 
structure akin to that of the cross-section of a flux tube that can 
be twisted into a spheromak in a 3D geometry (see MDD). The 
flux tube structure has a large inertia and is bound together by 
magnetic tension, possibly making it more resistant to dissipa- 
tion on kinetic scales than the smaller- scale current filaments. 
During the simulation the flux tube continuously gains mass and 
magnetic field until its further growth is limited by the periodic 
boundary conditions along the shock boundary. 

Comparing our work with the earlier ID simulations 
by |Bessho & Qhsawa| dl 999b and t he ID and 2D simula- 
tions (Dieckma nn et al.||2008| [2010), we find that the near- 
equipartiti on energy acce l eration predicted by these authors is 
confirmed. Mar tins et al.| ( |2009| ) in their 2D piston simulations 
found that the upstream ions ahead of the shock were filamented, 
extending the region of magnetic field growth. We find that the 
magnetic field growth is exponential in the foreshock region. 

Dieckma nn et al.| ( |2010| ) showed that the structures for a 
lower speed simulation and smaller magnetic-field-to-shock- 
normal angle are planar and ID. The simulations presented here 
show a greater departure from one-dimensional behaviour. The 
greater field angle and flow speed in our work allows more mo- 
tion transverse to the shock plane, which apparently triggers dif- 
ferent processes in the shock transition layer. 

The results shown here carry several important implications. 
Firstly, the filament generated magnetic field can be stored in 
magnetic vortices. We consider here a magnetic field that is rel- 
atively strong in that it yields an electron cyclotron frequency 
that is comparable to the electron plasma frequency. Secondly, 
oblique shocks have good acceleration properties, increasing 
electrons to near equipartition with ions. Thirdly, primordial 
fields are becoming increasingly accepted in GRBs and will have 
a decisive affect on the dynamics of plasma internal shock, ac- 
celeration of electrons and their injection int o the Fermi mech- 
anism. Fourthly, as is becoming evident (Lemoine & Pelletier 
|2010[ |Bret & Die ckmann 2010) a reduced mass ratio has an 
important effect on the results of PIC simulations and this ap- 
proximation needs to be validated by linear theory if used. More 
specifically in the context considered here, heavier ions can un- 
dermine the magnetic suppression of the filamentation instability 
that may work for ions with a low mass. 

Radiative los ses can be non-negligible for shocked colli- 
sionless plasmas (Fleishman & Toptygin 2007 ; Schlickeiser & 
|Lerche|2007||2008| ). Although the PIC framework does not take 
into account radiative processes, we can infer that at least the ra- 
diative cooling by synchrotron-type emissions should not affect 
the conclusions propounded in this paper. The ions do not ra- 
diate significant energy and only a fraction of the electrons will 
convert some energy to low-energy photons through synchrotron 
cooling, in particular because the individual particles only spend 
a short amount of time in the high magnetic field region. The in- 
teraction between the relativistic electrons and the photon seed 
that is not captured by PIC simulations may, however, result in 
stronger energy losses. 
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